Interaction of linear waves and moving Josephson vortex lattices in layered 

superconductors 



O 

(N 



3 



A. V. Chigine\0 and V. V. Kurin 

Institute for Physics of Microstructure of the Russian Academy of Science, GSP-105, 603950 Nizhny Novgorod, Russia 

(Dated: January 3, 2010) 

A general phenomenological theory describing dynamics of Josephson vortices coupled to wide 
class of linear waves in layered high-T c superconductors is developed. The theory is based on 
hydrodynamic long wave approximation and describes interaction of vortices with electromagnetic, 
electronic and phonon degrees of freedom on an equal footing. In the limiting cases the proposed 
theory degenerates to simple models considered earlier. In the framework of the suggested model we 
undertook the numerical simulation of resistive state in layered superconductors placed in external 
magnetic field and demonstrate excitation of linear waves of various origin by a moving vortex 
£^ , lattice, manifesting in existence of resonances on current-voltage characteristics. 

I—* 

CO . I. INTRODUCTION 

High-T c superconductors (HTSCs) with strong anisotropy has been the subject of intensive investigation during 
many years. The layered structure, intrinsic Josephson effect, and complex chemical composition provide a variety 
of physical properties of such materials. A great deal of attention has been paid to Josephson dynamics of layered 
5h ' superconductors. 

The application of an external magnetic field parallel to the layers to a layered HTSC leads to formation of a 
Josephson vortex lattice (JVL) which can be moved under the influence of an external current applied perpendicular 
to the layers. When this lattice collides with the edge of the structure the electromagnetic radiation is generated. 
' This principle may be used for construction of small and efficient electromagnetic oscillators with the frequency 
being limited from above by the energy gap in HTSC which is of the order of 10 THz. Recently, the radiation with 
the frequency of 0.85 THz and power of 0.5 /iW has been obtained from the structure based on F^S^CaC^Os 1 . 
The small and compact sources of electromagnetic radiation of terahertz band are very interesting for applications 
in astrophysics, medicine, biology and many other branches of science. Therefore, the investigation of Josephson 
dynamics in layered HTSCs is highly important. 

Layered HTSCs represent complex physical system with wide spectrum of eigenwaves, which includes electromag- 
netic, plasma, different phonon modes and Carlson-Goldman mode^. All these modes are coupled to vortices and can 
be excited by the moving JVL. The reverse action of the excited waves on JVL resulting in changes of vortex shape 
■ and mutual arrangement of vortices will affect the current-voltage characteristics (CVCs) of layered superconductor 
and intensity of electromagnetic radiation. Thus coupling of linear modes and moving vortices can be practically used 
(y) . for controlling the dispersion properties of eigenwaves by the external magnetic field, for diagnostics of eigenwaves by 
methods of Josephson spectroscopy^, and for controlling electromagnetic radiation from layered superconductor. 

The model adequately describing vortex dynamics in layered HTSCs should take into account the interaction of 
vortices with all weakly damping linear modes. However, the development of the theory of such comprehensive kind is 
not a simple task and a main obstacle is a nonlocality of the constitutive equations of a layered superconductor or, in 
other words, spatial dispersion. Even in hydrodynamic approximation when the dynamical equations are differential, 
the consistent approach to description of layered media requires solution of differential equations for field and matter 
in the domains of homogeneity, and subsequent joining of the solutions found for adjacent domains on interfaces. 
If the number of modes in the system is rather high then the exact description leads to sophisticated problem for 
eigenvectors and eigenvalues of high-dimension transfer matrix. Example of calculations of this kind applied to normal 
nonsuperconducting layered media can be found in Ref. 0- The description of layered systems can be considerably 
simplified by employing the long-wavelength limit which deals with smoothed dynamical variables averaged over the 
spatial period of the layered structure. The goal of the present paper is the development of such a theory for a layered 
superconductor with intrinsic Josephson effect. We propose an averaged, sufficiently simple, hydrodynamic theory 
accounting for a wide spectrum of linear waves and vortex degrees of freedom. The model is easily extendable for 
including additional linear modes and mechanisms of their exitation. 

By present, a number of models describing Josephson dynamics of layered HTSCs in some special cases has been 
proposed. The first and most known one is the local model accounting only for the magnetic (inductive) coupling 
between adjacent junctions of the stack which was applied to description of the artificial multilayer Josephson struc- 
tures 5 and later to layered superconductors with intrinsic Josephson effect' 1 . The charge coupling between the adjacent 
junctions in layered superconductors has first been investigated in Ref. for the case of spatially uniform distributions 
of Josephson phase. There have been also some attempts to combine magnetic and charge couplings into one model^. 
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The "global" coupling of the junctions via the external waveguide connected parallel to the long Josephson junction 
stack has been considered in our previous paper—. The influence of nonequilibrium effects, such as quasiparticle 
imbalance, to the Josephson dynamics of intrinsic junctions in layered HTSCs, has been investigated in Refs. Illl4l4] . 
The effects of the in-plane dissipation to the Josephson vortex motion in layered HTSCs has been studied in Refs. Ifa - 
H7I . The role of phonons in Josephon dynamics has been considered in Refs. [Ull^ for the case of direct excitation of 
phonons by the electric field, in Refs. I20II21I for phonon excitation due to phonon-assisted tunneling. All these models 
may be unified into one theory describing the linear waves and Josephson vortices in layered HTSCs. 

In the present paper we formulate the comprehensive theory accounting for vortex interaction with linear waves of 
different physical nature which unifies the models mentioned above. The prescription to design a theory of such kind 
is the following. Instead of exact consideration of high-T c superconductor as a layered medium we consider it as a 
continuous anisotropic medium which consists of superconducting and normal electrons and of several sort of ions. To 
describe the dynamics of such a media coupled to electromagnetic field we use Maxwell equations jointly with two- 
liquid anisotropic hydrodynamic equations for electrons and equations of adiabatic Born-Oppenheimer approximation 
for ions. But this set of equations describe only linear properties of layered superconductor in long-wave limit. To 
introduce vortex degrees of freedom to the model and allow finite jump of Josephson phase difference on dielectric 
layers we " discretize" the set of equations simultaneously restoring the nonlinear term with sinusoidal dependence on 
Josephson phase difference in interlayer current. In various particular cases our model is reduced to the known models 
considered earlier. In the framework of the suggested model we undertook the numerical simulation of resistive state 
in layered superconductors placed in external magnetic field and demonstrate excitation of linear waves of various 
origin by a moving vortex lattice, manifesting in existence of resonances on current-voltage characteristics. 

The paper is organized as follows. The next section is devoted to formulation of effective averaged model describing 
linear waves and vortex degrees of freedom in layered superconductors with intrinsic Josephson effect and to estimation 
of effective parameters of the model. The third section is devoted to numerical experiment showing the effects of 
excitation of linear waves by Josephson vortex lattice moving in a layered HTSC. In the conclusion we summarize the 
main results obtained in the paper. 



II. THE PHENOMENOLOGICAL DESCRIPTION OF JOSEPHSON DYNAMICS OF A LAYERED 

HTSC — HYDRODYNAMIC APPROACH. 



In this section we derive the set of equations describing Josephson dynamics of a layered superconductor treating it 
as a continuous but anisotropic medium. First, we write down the Maxwell equations and anisotropic hydrodynamic 
equations describing superconducting condensate, electronic quasiparticles and phonon degrees of freedom of super- 
conductor. As a result of this stage we obtain a system describing linear waves in a layered HTSC in a long wavelength 
limit. Then, substituting the derivatives over the coordinate perpendicular to the layers by finite differences and, si- 
multaneously, restoring the nonlinear expression for the interlayer supercurrent in Josephson form j s = j c sin 9 n , we 
obtain the desired set of equations. 




X 



FIG. 1: Side view of the multilayer Josephson structure. The numeration of superconducting and insulator layers, and the 
coordinate system, are shown. 

In this paper we adopt the widely used concept that high-T c superconductor can be treated as a sequence of 
superconducting layers coupled through dielectric layers by tunnel effect providing the existence of Josephson effect 
and quasiparticle conductivity in c-axis direction. Such a layered superconductor with chosen coordinate system is 
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shown in Fig. [T] The x axis is chosen laying in ab plane of a superconductor which is assumed to be isotropic in this 
plane, direction of the y axis is chosen coinciding with the direction of the magnetic field which is assumed to have 
only one component, both static and alternating, and the z axis is directed perpendicular to the layers along c axis of 
layered superconductor. For such magnetic field the electric field and other physical variables would not depend on y 
and the problem becomes two-dimensional one. The superconducting and adjacent dielectric layers are enumerated 
by number n, as shown in Fig. [TJ 

We start the derivation from the assumption that, first, the characteristic scale of the solutions is much greater than 
the layer structure period, and, second, the module of the Josephson phase difference over dielectric layer is much 
smaller than unity. These assumptions allow us to treat layered superconductor as a continuous medium, and, thus, 
to formulate the desired model in the continual limit. Then, after " discretizing" the continuous model, we will restore 
possibility of finite jump of the phase of the order parameter and finally obtain the set of finite-difference equations 
representing the desired model allowing vortex solutions. 

As the basic equations we use Maxwell equations for electromagnetic field together with hydrodynamic equations 
describing contributions from superconducting electrons, normal electrons, and phonons. The next subsections are 
devoted to derivation of constitutive equations needed to construct the desired model. 



A. Field equations and superconducting electrons 

In this subsection we derive the set of equations phenomenologically describing Josephson dynamics of a layered 
superconductor in a continual limit. To do this, first we use Maxwell equations together with anisotropic hydrodynamic 
equations for supercondicting electrons. Aftewards we supplement the set of obtained equations with constitutive 
equations describing contributions from normal electrons and phonons, which are to be derived in the next subsections. 

As the starting point of the derivation of the equations describing the superconducting electron subsystem we use 
the expression of the superconducting momentum which we denote as p via the phase \ of superconducting order 
parameter 

p = hV X +-A, (1) 
c 

here A is the vector potential of the electromagnetic field and the negative electron charge is explicitly taken into 
account. On this stage we suppose that there are no vortices in the material, so that [V x V%] = 0. Accounting for 
this, we obtain 

B=^[Vxp]. (2) 
Substituting this expression into Maxwell equation [V x B] = 47rc~ 1 j fot + c _1 E we get the evolutional equation for E 

1 • c 4-7T 

-E= — VxVxp jtot- (3) 

c 2e c 

The total current j tot in is the sum of the supercurrent j s , current of normal electrons j„, and ion current jj. 
The supercurrent is related to velocity by usual expression j s = — enov but in anisotropic Ginzburg-Landau model 
the momentum is connected to velocity by relation p = 2m s v via tensor m s of effective mass. We write down this 
tensor in the form m s — diag(m/r s , mT s ), where m = [ra xx m zz ) l l' 1 is the geometric average of the effective masses 
of superconducting electrons in the in-plane and interlayer directions, factor T s — {vn xx jm zz ) l l' 1 is a measure of 
anisotropy of superconducting properties of the material. Using the expression for supercurrent, one rewrites in 
components 

1. c fd 2 Pz d 2 Px \ An ( en T s \ 

1^ c f d 2 p z | d 2 p x \ 4tt / en 
dx 2 dxdz ) c \ 2mT 6 

Though we assume that nothing depends on y-coordinate in our system, in the Eq. (I4bp we take into account that 
the z-component of [V x B] actually consists of two terms 



c Ez = 2e + ~dx~d~z) ~ ~ [~ *rf?> + jn ' + *« " ^ ' ' (4b) 
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with the last term in the r. h. s. of this equation being proportional to the bias current j ex t- The Eqs. (|4a|) . (|4b[) are 
to be supplemented with evolutional equations for p. 

If we differentiate Eq. (JTJ) with respect to time and introduce chemical potential of superconducting electrons 
2/i = h\ — 2etp then we come to equation of hydrodynamic type 

mv + V^i = -eE. (5) 

Further we make a model assumption that chemical potential of superconducting electrons is defined by formula 
of degenerated Fermi gas fx = Bp where £p — h 2 (3n 2 n) 2 / 3 /2m is the Fermi energy and n is a concentration of 
superconducting electrons. Assuming that deviation of the concentration n s from its equilibrium value no is small we 
can linearize equation ([5]) what yields 

2 n - 

m s v s + -£f, s — = -eE, (6) 

3 n 

where we denote as £p <3 — ep(n — n ) a Fermi energy for unperturbed superconducting concentration. Further, 
excluding the concentration deviation n s from hydrodynamic equations using Maxwell equation which we write in the 
form (V • D) = Anp s = — 47ren s , we come to evolutional equation for p 

1 p = E-r2v(V-D), (7) 



or, in components 



2e 



1 . „ 2 fd 2 D x d 2 D z \ 

1 • p ^ (d 2 D d 2 D z \ 

Here = (2/3) 1 ^ 2 vp/uj p is a screening length of longitudinal electric field, v 2 F = 2ep/m, and uj 2 — 4ne 2 nom~ 1 are 
averaged Fermi velocity and plasma frequency respectively which are defined as vp = (y 2 Fx v 2 F z ) x l 2 , u> 2 — (uj 2 lo 2 ) x / 2 . 
Vector of electric displacement field D is defined so that it includes charges of normal electrons and ions as bound 
charges, i.e. 

D = E + 4vrP„ + AttV, = (/ + 47rx„ + 47TXi)E, (9) 

where P n and are the polarizations associated with normal electrons and phonons, respectively, so that 
(V • P n ,i) = —p n ,i- Tensors X n , \i are normal electrons and phonons susceptibility respectively. 

A set of equations Eqs. (|4"aj) , (|4b| . (|5a|). (|8b|) describes linear dynamics of electromagnetic field and superconducting 
condensate in layered superconductor in the long wave limit. Other degrees of freedom as normal quasiparticles and 
phonons of different kind come in the Eqs (|4a|). (|4b|) via extra current densities j„, ji and in the Eqs. (|5a|) . (18b[) 
via diplacement vector D. In the next two subsection we find contributions from quasiparticles and phonons to 
susceptibility and make the set of equations complete. 



B. Contribution from normal electrons 

To describe normal electrons we assume that they form degenerate Fermi gas and apply hydrodynamical Thomas- 
Fermi approach what gives the following equations 

n n + n 0n (V ■ v„) = 0, 

2 ef 

m n {-v n + vw n ) + Vji„ = eE. (10) 

3 n 0n 

In the Eqs. (fT0|) rh n = diag(m/r„, mT n ) is the mass tensor of normal electrons, T n is the anisotropy parameter, no n 
is the unperturbed electron concentration, n n is the deviation of electron concentration from its unperturbed value, 
v„ is the electron velocity, ep „ = ep(n = no.„) is the Fermi energy of normal electrons and tensor v — diag(i^;, v z ) 
characterizes the electron collision frequencies. For superconductor the condition u) <^ v x _ z is usually satisfied and 
equations can be somewhat simplified by neglecting the term v in comparison with vx. Then such hydrodynamic 
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equations can be solved and the relation between normal current j n and electric field E can be easily obtained. In 
(u>, k)-representation one gets 



w 2 I iiujv z T n ~2v 2 Fn k 2 z 2v 2 Fr k x k z 
4ttA 



2v Fn k x k z 3iujv x /T n - 2v 2 Fn k 2 x 



E = -icJX«E. (11) 



Here w 2 n = 4ire 2 no n /m is the plasma frequency, vf u — ( £ f /m) 1 / 2 is the Fermi velocity, both related to normal 
electrons, A = Ziujv x v z — 2v x T~ l v Fn k 2 — 2v z Y n v 2 Fn k 2 . Without spatial dispersion of normal electrons, i. e. at 
vp n = 0, this equation became the simple anisotropic Ohm law 

j n = <tE, a = ding(a xx ,a zz ) 7 (12) 

where a xx = ui 2 n T n /(4iris x ), a zz = u> 2 n /(4Tru z T n ) are respectively the in-plane and interlayer normal conductivities. 
Here it is worth to mention that such local model of normal conductivity has been used in works — ~— devoted to 
study the influence of the in-plane normal conductivity to the dynamics and stability of moving JVLs in layered 
HTSCs. Our approach naturally introduces into consideration the spatial dispersion originated from nonzero pressure 
of normal electrons. 



C. Phonon contribution 



This subsection is devoted to accounting for phonons in our phenomenological model describing Josephson dynamics 
of layered HTSCs. In the present paper we take into account only the direct excitation of infrared-active phonons by 
the electric field of Josephson oscillations. In fact here we will actually generalize the approach used in Refs. [lafl9l 
to the case of non-uniform solutions in distributed Josephson systems. The nonlinear interaction between Josephson 
oscillations and phonons due to effects of phonon assisted tunneling considered in Refs. I20II21I would not be considered 
in this subsection, but later, when we formulate general nonlinear equation, we will show how this mechanism can be 
introduced in our model. 

In order to find the phonon contribution to the dielectric permittivity we use standard approach (see, for example, 
ReflU). Let us write the classical equation of motion for ions derived in Born-Oppenheimer adiabatic approximation, 

M„z£ = - <%£ M z& + q„E^. (13) 

In this formula z^ is the ion displacement from the equilibrium position, v is the ion index in the unit cell, N = 
ni&i + n 2 a 2 + n 3 a 3 is the unit cell index, a 1:2 ,3 are lattice periods, G^_ m is the "seed" force tensor resulting from 
the interaction between ions via valence electrons, q v is the ion charge, M v is the ion mass, Ej^j is the microscopic 
electric field at the point of the ion. The interaction between ions and conductivity electrons is carried out via the 
last term in the r. h. s. of (fT3|) . In the equation (fl~3|) and further the product of a tensor and a vector is written in 
components as (Ax.)i — AijXj. The relation between microscopic field and average macroscopic field can be written 
via so called Lorentz tensor L^_ m 

E£j = E N + 4tt Y, ^n-mPm, (14) 

where Pj^ = g^z^j is the dipole momentum of the /i-th ion in the M-th cell. With the account for the difference 
between microscopic electric field from macroscopic field the equation for ^-th ion motion takes the form 

M V Z£ =-J2 ^N-M Z M + ^ E N, (15) 

where F-^^^ = — / i'nq u L^ i m q ti is renormalized force tensor. 

The contribution of ions into dielectric permittivity s may be found from the expression for ion current density 



jN = ^XX i; N> (16) 



6 



where V is the unit cell volume. After some algebra we get the ion current density 



V 2. _,„2+,„2. rif ^ V M.a*».. = - lUXi ^ (17) 



where a is the phonon mode index, L is the number of ions in the unit cell, w p ^(k, a) is the phonon frequency, ei,(k, a) 
is the polarization vector. w p ^(k, a) and e„(k, a) are yielded from the equation for eigenvalues and eigenvectors for 
the matrix F v "(ls.) 

^F^(k)e M (k,a) = M v GJ 2 ph (k,a)e v (k,a), (18) 

A" 

where F^k) is the Fourier image of the tensor F^^ M 

F^(k) = e" lkN F^- (19) 

N 

The expression (I17p is the well-known formula describing phonon contribution to dielectric permittivity, which may 
be found in many solid state physics courses. Of course, the concrete definition of phonon frequencies and polarization 
vectors for layered superconductors require detailed spectroscopic or numerical investigation. Example of calculations 
of phonon properties of typical layered superconductor with intrinsic Josephson effect Bi 2 Sr 2 CaCu2 08 can be found 
inRef.[H. 

So now we have found the contributions from normal electrons and phonons to the dielectric permittivities and 
are able to define electric displacement vector D = (I + AitXn + 47TXi)E = eE, using expressions (fTT|) and (|T7|) for 
susceptibilities \ n and Xi- 

The Eqs. (ga]), (|4b|). (|5a|) . ([8b| together with the relations (fTTj) and ([TT]) represent complete system describing, in 
principle, all linear waves in long wavelength limit. However, this system does not describe Josephson vortex degree of 
freedom in layered HTSCs because it is linear and does not allow the finite jumps of the phase of order parameter. In 
the next subsection we will show how Josephson vortices can be included in the model and formulate general nonlinear 
model describing both Josephson vortex dynamics and linear waves in layered HTSCs. We will compare this model 
with several known ones, and reveal the relations between the parameters of continual model and the parameters of 
layered superconductors. 



D. Discretization of the model 



The aim of the present subsection is the transformation of the set of equations in the continual limit, to the form 
admitting solutions in the form of vortices, which are able to move along the layers of the HTSC. To do this, we 
split the continious medium into the series of layers of thickness s equal to the period of the layers of the HTSC, 
and oriented perpendicular to the c axis of superconductor. We introduce new variables describing the corresponding 
fields at the moment t, coordinate x, and in some point within the n-th layer. It is possible provided field distributions 
are smooth inside the layer. After such discretization we allow the phase of order parameter to make a finite jump 
between layers and restore the nonlinear expression for the interlayer Josephson current. 

First of all, let us decompose all dynamical variables to two sets defined on two different lattices shifted with respect 
to each other. The nodes of these two lattices can be thought to be placed in the middle of superconducting and 
dielectric layers, respectively. The nodes of these sublattices we will denote by the number of period n. We attribute 
the variables E x , D x ,j x ' s ' l ,v x ' s ' 1 to the first set and E z ,D Z} j™' s ' l ,v™-' s ' 1 to the second one. After such a decomposition 
we replace the derivatives of the dynamical variables over the coordinate perpendicular to the layers in the set of 
Eqs. flU}, (gb|, (JHaJ) , (j8b}, (HI]), (17]), by finite differences, using the following rule : 

dV ^d 2 V 
s— ^±V n±1 TV n , s 2 —~y -> V n+ i - 2V n + V n -i = A n V n , (20) 

where V n stands for the one of the dynamical variables of the problem, upper sign should be used for variables from 
the first set and lower sign — for the second one. Such a rule of discretisation follows from integral form of Maxwell 
and hydrodynamic equations and provides necessary symmetry of finite-difference equations. To demonstrate how this 
procedure works let us apply it to well known telegraph equations describing electromagnetic waves in transmission 
lines 



Li + U x = 0, CU + I x = 0, 



(21) 
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here U, I are voltage and current in the line, L, C are linear densities of inductance and capacity, respectively. 
Attributing voltage and current to different sets and replacing spatial derivatives by finite differences using our rule 
we come to equations 

sLi n + U n - U n - X = 0, sCU n + J n+ i - /„ = 0, (22) 

expressing two Kirchhoff laws for discrete L, C chain. Now let us return to our problem. Applying the formulated 
procedure to our set of differential the substitution (f2TJ|) we obtain the following system 

I^tti c d . . c Air f en Q T s , n .,- t \ 

cdt xn = 2^sdx {Pzn ~ Pzn - l} ~ 2e^ nPxn ~ ~c \ 2rrT Pxn+Jxn+3xn ) > (23) 
I d _ c d 2 c d Air ( en Q .„ . 4 . \ 

cdt zn ~ ~2edx* Pzn ~ 2e~sdx~ [Pxn+1 ~ Pxn) + ~ \2mF s Pzn + 3zn + 3zn ~ 3ext ) ' ( ' 

~hwt Pxn = Exn ~ r * {h Dxn + \ik {D * n ~ D * n - i} ) ' (25) 

15 ./IS, 1 \ 

~2edt Pzn = Ezn ~ rd {sdx {Dxn+1 ~ ° xn) + f 6 * **) ' (26) 

D = eE, j"-' = ^(Xn,iE), i = i + 4nXn + ^Xi- (27) 
at 

The Eq. (|27j) of this system combines the contributions from normal electrons and phonons which are given by the 
formulas (jTTj) and (fTT|) for corresponding tensors of susceptibilities. In a common case, when spacial dispersion 
takes place, these tensors contain wavevectors in the interlayer direction. Therefore, the procedure of discretization 
should be applied also to Eq. (|27p accounting for specific expressions for normal electron and phonon contributions 
to dielectric permittivity. 

Now this system allows finite jumps of physical variables but it still remains linear and, therefore, does not describe 
Josephson vortices. To describe them we need to introduce Josephson nonlinearity in the expression for interlayer 
supercurrent. This procedure is made in the following way. Integrating the expression for z-component of the 
superconducting momentum over z from n-th to n + 1-th layer, we get Josephson phase difference 

n = -- p z dz = Xn - Xn+i ~ / A z dz. (28) 

Then we need to substitute 6 n into the set of equations in a correct way. As the interlayer supercurrent in layered 
HTSCs has Josephson nature, we need to substitute p Zn — > — Hs -1 sin 6 n in the last term of the Eq. (|24[). The 
remaining terms in Eqs. (|23[) and (|24j) having p Zn are in fact the components of [V,B], that is why in these terms 
p Zn — i —hs~ 1 6 n . The Eq. (|2l)j) at rn — is actually the Josephson relation H9 = 2eU, then in this expression also 
p Zn — > —hs~ 1 8 n . At rd the Eq. (126)) describes violation of the Josephson relation due to charge effects; the 
influence of these effects on Josephson dynamics is considered later, in the section devoted to numerical experiment. 
Let us rewrite the system (f2"3")l , (j2"4")l . (|2"5|) . (I2U1) . (|27p in the notations used when describing bulk superconductors: 

$n ^p s 2 _ c$ _ c 2 2 v 2 Fs d„ 

a„k — 



c 8n 2 cs r s ' ab 8Tr 2 j c sT 2 T s cu 2 ' d ui 2 s 

" Jt - s " P s Ps 

As a result, we get the following set of equations: 

cdt Exn ' ~2^ 2 dx {6n ~ 9n - l] + T^ A " J - " ^ xn +3 Xn +J xn ), (29) 



cdt " 2irsdx 2 cs dx '• 



= " — "-jr-dln+l ~3in) ~ ~0c + " jext), (30) 



— gi^n - E xn - -j- [g^D xn + --(£>, B - D, n _ x ) j , (31) 

itsWt dn = Ezn ? - + i A -^«) • (32) 

D = eE, j^ = ^(^E), e = / + 47rx„+47rx l . (33) 



The model represented by these equations (I3TJ)) . (f3"Tj) , fl32}, (133)) describes interaction of Josephson vortices 

with electromagnetic waves, plasmons, and phonons on an equal footing. This model possesses the necessary symmetry 
with respect to permutation of x, z coordinates and describes the spatial dispersion resulted from electron and ion 
degree of freedom in the system. Besides, our system allows to account for the influence of linear modes of any nature 
to the dynamics of Josephson vortices in a layered HTSC, by inclusion of the corresponding terms into the dielectric 
permittivity e. Say, in recently discovered FeAs based superconductors 2 ^ - — it can turn out that magnetic degrees of 
freedom play an important role. Our approach allows to include them in this general scheme. 

Of course, our general model contain the particular theories considered in earlier works as a limiting cases. Among 
these special cases are the model with the magnetic coupling between the layers^, the model with the charge coupling 
between the layers^, the model taking into account the in-plane quasiparticle current, the model accounting for the 
infrared-active phonons polarized perpendicular to the layer a 18 ' 19 , etc. Results of two works 8 -^, where the attempt to 
unify magnetic and charge coupling have been undertaken also can be reproduced by our model. But our model have 
extra terms in Eqs. ([3"Tj) and containing d/dx which are absent in Refs§0. Although small far from resonances, 
these terms become significant when Josephson frequency get close to phonon and plasmon frequencies. 

The set of Eqs. (f!25|) . (13TJ1) . (f3~Tj) . (|3^|) . ([33)1 contains the only nonlinear term, which describes Josephson nonlinearity. 
While deriving this system we have neglected all hydrodynamic nonlinearities such as nonlinear dependence of a 
pressure on concentration, difference between Lagrange and Euler varyables, nonlinear expression for the current 
density and so on. The account for these terms would lead to the existence of combination processes, such as, for 
exaple, Raman and Brillouin scattering. Due to strong anisotropy of the layered HTSCs, the most significant effects 
of this type are the ones providing the dependence on concentration of the parameters determining the transport and 
elastic properties of the material in the direction of c-axis, such as Josephson critical current density j c , interlayer 
normal conductivity a zz , etc. The accounting for these dependencies leads to such effects as phonon-assisted tunneling 
due to the interaction of the electrons with Raman-active phonons. This effect has been considered in Ref. [2l| for 
acoustic phonons in long Josephson junction and later in Ref. |2Q. for intrinsic junctions. The similar effect should be 
also expected for plasmons. In the present work we restrict ourselves to the case when these effects are unimportant. 

Let us consider now the examples how some known models arise from our equations in the limiting cases. For 
instance, if we set the phonon current, the in-plane displacement and normal currents to be equal to zero {j x ,z ph — 
dE xn /dt = j™ n = 0), propose the interlayer normal current to have the purely Ohmic character (j Zn = a zz E zn ), and 
neglect the nonzero screening length of the longitudinal electric field (r<2 = 0), we obtain 



o 2 e n ( i . 1/ 2 d 2 e n , $ do n , . „ j ext 




8x 2 \X 2 " X 2 C \ 3 dt 2 zz 2ircs dt " j c 



(34) 



where Xj = T s s, X c = r s A a b, oj 2 — (87r 2 j c cs)/$o- First this set of equations has been derived in Refs.H@ and describes 
the Josephson dynamics of stacked distributed Josephson junctions and layered superconductors with the magnetic 
coupling between the layers. In another case, when we suppose the lateral dimensions of the system to be smaller 
than Josephson length, so that d/dx — 0, neglect phonons, and set j Zn — a zz E zn , we obtain 

1 d 2 6 n 1 / j ext \ <±ko zz d0 n 

= -t^(1-t?A„) sm6» n -— ~t 5j~ ' ( 35 ) 



c 2 dt 2 X 2 \ j c J c 2 dt 

where r\ = r^/(sd s ) is the parameter of charge coupling. This set of equations has been first considered in Ref. and 
describes the dynamics of the layered superconductor with the charge coupling between the layers. Other particular 
examples of models considered earlier may be obtained from our model in a similar way. 



III. DISPERSION CHARACTERISTICS OF LINEAR WAVES. 



The layered HTSCs has a wide spectrum of eigenwaves which, in principle, may be radiated by a moving JVL. The 
linear modes excited by JVL, in turn, may affect the shape and mutual arrangement of vortices in the lattice. The 
excitation of linear waves by the moving JVL leads to appearance of resonant steps on CVCs with the frequencies 
being equal to the ones of the radiated modes. In order to identify the resonances on CVC it is necessary to know 
the dispersion characteristics of linear waves in layered HTSCs, which can give much information about eigenmodes 
in the material and conditions of their excitation. 

In order to build the dispersion characteristics we first linearize the set of Eqs. (|29]). (J30J) , ([3"T]) . (j32]) . (f33j) assuming 
\6 n \ <C 1, so that sin6*„ w n , and set the bias current to zero. We also neglect dissipation in the system (a xx = 
a zz = 0) keeping in mind that it may be accounted for by perturbation theory. Taking all dynamic variables 
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~ exp (ikx + iqn — iuit) and setting the determinant of the obtained linear system to zero, we get the dispersion 
equation of the linear modes. It has the standard form 



det 



-£totaiij{u, k) + kikj - k 2 5i 



= 0, 



(36) 



where itotai contains linear contributions from superconducting and normal electrons, and phonons. For short, we do 
not write this equation in an explicit form. We also assume the simple model expression of phonon susceptibility 



Xph 



1 

4-7T 



V 









o 2 



'ph 



(37) 



u 0z 



This expression describes two optical phonon modes: one polarized along the layers (x-phonon) and other polarized 
perpendicular to the layers (z-phonon). Here loq x and loq z are the frequencies of the x-phonon and the z-phonon, 
respectively, 7^ are 7^ the damping coefficients of these modes, Q x and VL Z are the oscillator strengths. For 
simplicity we assume wo IjZ , lph x z , and Cl x>z to be independent of the quasimomentum, i. e. the bare phonon modes 
has no spatial dispersion. Moreover, in this section we neglect phonon damping (~/ph , = 0), assuming that it may 
be accounted for by perturbation theory. We will return to the nonzero jph _ later, in numerical experiment. 

In our previous work we have considered wave vector surfaces for the analysis of the dispersion characteristics of 
the layered superconductor with phonons in the continual limit^. Now we build the dispersion curves u>(k,q) for the 
periodic system. 

The dispersion curves of linear modes in a layered HTSC are schematically shown in Fig. [5J The electromagnetic, 
plasma, and phonon modes are shown. We build the dispersion characteristics for certain directions in the Brillouin 
zone. The curves in the right part of the plot are for the T — X direction (growing k), the curves in the central part of 
the plot are for the T — X direction (growing q) , and the curves in the left part are for the X — W direction (growing 
k), along the edge of the Brillouin zone. We use standard notations T, X, W for the points of the Brillouin zone. The 
dispersion curves for the case of absence of the interaction between electromagnetic and plasma modes with phonons, 
i. e. at zero oscillator strengths (Q x ,z = 0), are plotted by dashed lines. It is seen that at nonzero oscillator strengths 
the splitting of dispersion curves appears in the places of intersection between bare curves. 

We note that due to anisotropy the separation of modes into the electromagnetic and plasma mode is relative, 
depending on the direction of the wave. For example, consider the lowest branch of the dispersion characteristics. 
The wave having q = is the electromagnetic mode and the one having k = is the plasma mode. And vice versa, 
for highest branch, the wave with q = is the plasma mode and the one with k = is the electromagnetic mode. For 
the arbitrary direction of the wave vector is it impossible to say whether it is an electromagnetic or plasma wave. 

The frequency of the lowest branch of the dispersion characteristics corresponding to the point T {k = 0, q = 0) 
is ujj — the frequency of the Josephson plasma resonance. In the notations of the hydrodynamic model it is equal to 
w p s /vTV The typical values of Wj for layered HTSCs are of the order 100 GHz. At the values of the parameters 
corresponding to layered HTSCs, the frequencies of the highest branch of the dispersion characteristics are in the 
optical range. Therefore, this mode cannot affect the JVL motion in the layered superconductors. 

The important particular case of the modes in a layered HTSC is the Swihart wave, i. e. the mode belonging to 
the lowest branch of the dispersion characteristics, which propagating along the layers and having the standing wave 
structure in the direction perpendicular to the layers. The dispersion curve w(fc) of this wave has the form of the 
hyperbole. The slope of the asymptotes of this curve determines the Swihart velocity 



c 2 {q) = n- 



-} + 27^(1- cos q) 
1 + 2/i(l — cos q) 



(38) 



where q is the transverse wave number of the Swihart wave. It is seen that the antisymmetric (q = ir) Swihart mode 
is the slowest one, and the symmetric (q = 0) mode is the fastest one. We will use the formula (|38|) for the analyzing 
the results of the numerical experiment. 



IV. NUMERICAL EXPERIMENT 



Now we apply the derived system |29|) . (f30|) . (f3lj) . (f32j) . (|33|) to the numerical investigation of the dynamics of 
Josephson vortices in layered superconductors with the account for their interaction with various linear modes. The 
motion of Josephson vortex lattice (JVL) leads to excitation of linear modes of a layered HTSC which, generally, affect 
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FIG. 2: (Color online) The schematic drawing of the dispersion characteristics of linear waves in a layered HTSC. The dashed 
lines show the dispersion without the interaction with phonon modes. The inset shows the fragment of the Brillouin zone of 
the layered structure. 

the moving lattice, leading to distortions in vortex shape and changes in their mutual arrangement. The excited linear 
modes may lead to resonances on CVCs of layered HTSCs with moving JVL. In the present section we numerically 
investigate the excitation of linear waves by the moving JVL, and their influence on CVCs, accounting for possible 
changes in vortex shape and their rearrangement. 

In the numerical experiment we use periodic boundary conditions in both in-plane and interlayer directions for all 
variables. However, for Josephson phase difference 9 n the boundary conditions in x-direction is modified so that 

9 n (L)=9 n (0) + 2irR ni (39) 

where R n is the number of vortices trapped in the n-th junction of the stack, and L is the length of the system in 
the longitudinal direction. Instead of 9 n we introduce 9' n = 9 n — 2itR n x/ L satisfying the periodic boundary condition 
in the longitudinal direction. Accounting for this, we write the expression for the interlayer Josephson current 
as j c sin(9' n + 2irR n x/L). The choice of boundary conditions described above provides the simplicity of numerical 
solution of the system; however, when using such conditions one takes into account only volume effects, neglecting 
the influence of boundaries on Josephson vortex dynamics. 

In order to solve the system (|2T)1) , ([3TJ)) , (pTTj) , (|3"2"|) , (j33"|) numerically we transform it to the set of evolutional equations 
in ordinary derivatives in time. To do this, we use the exponential Fourier transform by standard 2D-FFT algorithm. 
The obtained set of ordinary time-dependent differential equations is solved by Krank-Nikolson scheme. The similar 
approach to the numerical experiment has been used in Ref. [Tol . 

The periodic boundary conditions along the layers imply that the number of vortices captured in each junction of 
the system is a constant. In our calculation we assume the number of vortices to be the same for all layers and denote 
it as R. 

In the numerical experiment we use the typical values of the parameters of Bi2Sr2CaCu20s: j c — 150 A/cm 2 , 
\ a b = 1700 A, s = 15 A, a a i, = 5 • 10 4 (ft- cm) -1 , a c = 2 • 10~ 3 (f2- cm) -1 , e zz = 12. The parameters of phonons have 
been taken from the experimental data obtained by spectral ellipsometry 2 ^. In the following we assume the length of 
the system L to be measured in the units of Josephson length Xj — T s s and the bias current density to be measured 
in the units of Josephson critical current j c . In the calculations we use values of L corresponding to the ones of the 
patterns used in experiments. The number of vortices in a layer is set so that it corresponds to the dense JVL (at a 

— 1/2 

given L). We also assume the frequencies and the voltages to be measured in the units of uj p T s 

A. Results 

This subsection is devoted to the results of the numerical modeling of the dynamics of the JVL in a layered HTSC 
performed basing on the derived model. First, we investigate the simplest case, when the spatial dispersion, phonons, 
and all interlayer couplings except the magnetic one are neglected. Then, we complicate the problem step-by-step, 
introducing phonon susceptibility, charge coupling, and study the effects caused by these complications. 
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1. The simplest case — the magnetic coupling without phonons and spatial dispersion. 

To check the efficiency of the used numerical scheme we perform the calculation of the CVC of the layered HTSC 
with the moving JVL in the case of absence of spatial dispersion, phonons, and in-plane dissipation. The interlayer 
normal current is assumed to be purely ohmic. This approximation corresponds to the case of magnetic coupling 
between the layers, which has been considered in Refs. yfl 

Consider the CVC built for the case R — 6 which corresponds to typical values of the external magnetic field 
(Fig. [3]). At j ext — a triangular JVL is established in the structure (Fig. 2]). As the bias current increases, the 
JVL moves with the increasing velocity, keeping the triangular vortex arrangement. When the bias current reaches 
0.18, the JVL velocity stops growing, so that the voltage on the structure is established at the value 0.18. This step 
denoted by 1 appears due to coincidence of the JVL velocity with the characteristic velocity of the antisymmetric 
(q = 7r) Swihart mode of a layered HTSC. 

The JVL remains triangular on the whole step 1; the vortex shape exhibit specific distortions due to the resonant 
increase of the harmonics with q = ir which grow as the bias current increases. The distribution of the magnetic field 
in this regime does not qualitatively differ from the static one (Fig. |4j . 
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FIG. 3: A CVC of the layered HTSC with the moving JVL. The digits enumerate the steps of the CVC in the order of resonance 
frequencies increase. The insets schematically show the mutual arangement of vortices on different resonance steps. 
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FIG. 4: (Color online) The magnetic field distribution corresponding to the static vortex lattice (j ex t = 0). For the sake of 
convenience, here and in the subsequent figures we add constants to the distributions of the magnetic field in different junctions. 
The actual value of the constant component of the magnetic field is as for the junction 1. 

At the value j ex t = 0.44 the CVC jumps from the step 1 to the Ohmic branch. The amplitude of the electromagnetic 
field on the Ohmic branch is small and the vortices form rectangular lattice. 

When the bias current reaches the value j ext ~ 1.0 the second step on the CVC appears. This step is due to the 
resonance with the symmetric (q = 0) Swihart mode. In this regime, the moving vortices still form the rectangular 
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FIG. 5: (Color online) The magnetic field distribution corresponding to the end of the step 2 in Fig. [3] Here j ext = 1.4. Giant 
amplitude of the electromagnetic field. 

lattice, but the amplitude of the electromagnetic field sharply increases. This regime of vortex motion is characterized 
by large amplitude of the electromagnetic wave which accompanies the moving JVL, the amplitude of the wave grows 
with the bias current increase (Fig. [5]). 

With further bias current increase the CVC again jumps to the Ohmic branch. The magnetic field distribution in 
this regime does not qualitatively differ from the one between the steps 1 and 2. 

The CVC obtained in our calculations is similar to the one obtained for two-stacked long Josephson junctions^. 
The only difference is that in our CVC the voltages of the steps differ by two orders from each other. This is due to 
the fact that the magnetic coupling in layered HTSCs is usually much stronger than the one of artificial multilayer 
structures. 

Consider now the CVC of the layered HTSC with the moving JVL at smaller external magnetic field (R = 4) 
(Fig. One can see that, in addition to steps 1 and 3 corresponding to the resonances with the antisymmetric 
Swihart mode (u = 0.14) and the symmetric mode (u = 32), respectively, there is a step 2 corresponding to half a 
frequency of the resonance with the symmetric Swihart mode. The magnetic field distribution on this step is shown 
in Fig. [7J It is seen that there are two oscillations of the magnetic field per spatial period of the system, though four 
magnetic flux quanta are captured in each junction of the stack. This regime also shows a large amplitude of the 
electromagnetic field. The distributions of the magnetic field and the mutual arrangement of vortices at the steps 1 
and 3 and between the resonances are similar to the ones for the case considered earlier (see Figs. |3l |4j [5]). 
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FIG. 6: CVC of the layered HTSC with moving JVL for weaker external magnetic field. 
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FIG. 7: (Color online) The magnetic field distribution corresponding to the end of the step 2 in Fig. [6] Here j ex t = 0.58. 



2. Excitation of a phonon by moving vortex lattice. 



The complex chemical composition of layered HTSCs provide large amount of phonon modes in such materials. 
Among these modes, there are "soft" phonons having frequencies of the order of several THz, which are smaller than 
the frequency of the energy gap in HTSCs. This makes possible the excitation of such phonon modes by a JVL moving 
in a layered HTSC. In this subsection we investigate the phonon excitation by a moving JVL. 

For the calculations we slightly complicate the model used in the previous subsection, introducing the sim plif ied 
expression ([57)1 for the phonon susceptibility. The phonon parameters used in calculations are taken from Ref. l28l . 
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FIG. 8: (Color online) The amplitudes of alternate components of the electric field and the polarization vs. external current in 
the vicinity of the phonon frequency. 

The Figs. [8] and [9] illustrate the excitation of the phonon z-phonon by the moving JVL. The Fig. [8] shows the 
dependence of complex amplitude modules of P z and E z harmonics with q = and k = kiattice, on the bias current. 
Here P z is the z-component of the phonon polarization and kiattice = 2ttRL~ 1 is the wavenumber of the main harmonic 
of the JVL. The peak on the dependence P z {jext) is due to the excitation of the phonon mode and appears at the 
value of j ex t which corresponds to the JVL motion with the frequency of the phonon. The Fig. |9] shows the peak on 
the contribution to the CVC due to the extra energy needed to excite the phonon. The similar phonon peaks on the 
CVC has been obtained earlie r 18 ' 19 for the case of spatially uniform Josephson junctions and junction chains. Our 
model actually generalizes the one used in Ref. to the case of distributed Josephson juctions and junction stacks. 

In the present simulations we do not consider the excitation of the cc-phonons. The reason is that the JVL is 
rectangular at low enough external magnetic field and at the frequencies close to the phonon one. As one can see from 
the Eqs. ([%!?)) , ([317)) . (j3"T)) . (f3"2|) . x-phonon is not excited by a rectangular JVL, as E x — in such a lattice. However, 
if the tensor of phonon susceptibility contains non-diagonal components, the excitation of x-phonons by E z of the 
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FIG. 9: Contribution to CVC due to excitation of the phonon mode in a HTSC. 

moving rectangular JVL is possible. 

In this subsection we have considered the excitation of a phonon by a moving JVL provided the condition that the 
phonon frequency is far from the frequencies of Swihart modes at a given ki a uice- The next subsection is devoted to 
the situation when the phonon frequency coincides with the frequency of one of the Swihart modes. 

3. Excitation of the hybrid phonon+Swihart mode by a moving Josephson vortex lattice. 

By choosing the external magnetic field applied to the structure, it is possible to make the frequency of the Swihart 
mode at k — ki attice equal to the phonon frequency. The Fig. [TT] shows the contribution to the CVC in the vicinity 
of the bias current value corresponding to the phonon frequency. Two peaks located close to each other and having 
nearly equal height appear due to excitation of two modes with close frequencies. To explain this effect, consider the 
fragment of the dispersion characteristic of linear waves in a layered HTSC in the vicinity of the point where the 
Swihart mode and the phonon mode interact (Fig. 1 10[) . In the absence of the interaction the dispersion characteristic 
would have the shape shown by dotted line, here the slanted line shows the dispersion of the Swihart mode and the 
horizontal line shows the phonon dispersion. 

In the presence of the interaction between two modes (Sl x , Vl z ^ 0) the dispersion characteristic takes the form 
shown by solid line in the Fig. 1101 The magnitude of the dispersion curve splitting is determined by the oscillator 
strength of the phonon mode (Q x or Q z , depending on the polarization of the mode interacting with the Swihart 
wave). The intersection points of the vertical line corresponding to ki at uce, and dispersion curves, give the resonance 
frequencies. 

The Fig. [TTj shows the contribution to the CVC from the excited hybrid phonon+Swihart modes. The distance 
between two peaks is equal to the difference between the resonance frequencies obtained from the dispersion charac- 
teristic. We note that the height of the peaks in Fig. [TTJ is of two orders higher that the height of the peak caused by 
the excitation of the pure phonon mode (Fig. [9]) . 

4- The violation of the Josephson relation — separation of normal electron charges by a moving JVL. 

In this subsection we investigate the influence of the charge effects on the JVL motion in layered HTSCs. Starting 
again from the simple model with the magnetic coupling, we now assume the parameter of the charge coupling to be 
nonzero r; 7^ and find out the differences in JVL dynamics compared to the case r\ = considered above. 

Consider CVC of the layered superconductor, calculated in the absence of phonons and in the presence of the charge 
coupling (Fig. [TJl dashed line). As it is seen from this figure, there are two steps of this CVC. As in the absence of the 
charge coupling, they appear due to resonance with the antisymmetric and symmetric Swihart modes. Let us consider 
the differences between CVCs in Figs. 151 and IT21 The first one is that the frequency of the first step is shifted towards 
higher frequencies, while the frequency of the second step is not changed. As in the case of zero charge coupling, the 
step positions correspond to the formula uj res — c(q)2TrRL^ 1 , where c(q) is determined by the expression (1551) . 

The second feature of the regime with nonzero charge coupling is that the amplitude of the second step is much 
smaller than in the case of zero charge coupling. To explain this effect, consider the dependence of the Josephson 
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FIG. 10: (Color online) Dispersion characteristic of the symmetric Swihart mode and phonon mode in the vicinity of their 
interaction point. 
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FIG. 11: Contribution to CVC due to the excitation of the hybrid phonon+Swihart modes. 

phase growth rate on the bias current (Fig. [T^J solid lines). It is seen that the phase growth rates vary from one 
junction to another or, the same, the vortex chains in different junction have different velocities. At the same time, 
the voltages on each junction are the same. According to the Eq. ([32]) . this is the demonstration of the violation of 
the Josephson relation. As it is also seen from Fig. [T^J only two or three junctions of possible four ones are locked to 
the second resonance of CVC. Therefore, the range of bias currents where the system remains on the resonant step, 
is smaller than in the case of zero charge coupling, when all junctions are locked. 

V. CONCLUSION 

We propose the comprehensive phenomenological model which describes the dynamics of the non-uniform distri- 
butions of Josephson phase difference in layered HTSCs, e. g. moving Josephson vortices and linear waves of any 
nature. Basing on this system we numerically build CVCs of a layered superconductor with the moving JVL and 
demonstrate the excitation of linear modes by moving vortices. The proposed model is shown to cover many effects 
which have been studied in previous works; in addition, we observe some new effects such as excitation of a phonon 
and hybrid modes by a moving JVL in layered superconductors. 
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FIG. 12: (Color online) Dependence of phase growth rate in Josephson junctions of the structure, on the bias current. Dashed 
line shows the constant component of the electric field E z in each junction. 
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